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ABSTRACT 

We study a large sample of 625 low-redshift brightest cluster galaxies (BCGs) and 
link their morphologies to their structural properties. We derive visual morphologies 
and find that ^ 57% of the BCGs are cD galaxies, ^ 13% are ellipticals, and ^ 21% 
belong to the intermediate classes mostly between E and cD. There is a continu¬ 
ous distribution in the properties of the BCG’s envelopes, ranging from undetected (E 
class) to clearly detected (cD class), with intermediate classes (E/cD and cD/E) show¬ 
ing the increasing degrees of the envelope presence. A minority 7%) of BCGs have 
disk morphologies, with spirals and SOs in similar proportions, and the rest 2%) are 
mergers. After carefully fitting the galaxies light distributions by using one-component 
(Sersic) and two-component (Sersic-I-Exponential) models, we find a clear link between 
the BCGs morphologies and their structures and conclude that a combination of the 
best-fit parameters derived from the fits can be used to separate cD galaxies from 
non-cD BCGs. In particular, cDs and non-cDs show very different distributions in the 
Rc-RFF plane, where Re is the effective radius and RFF (the residual flux fraction) 
measures the proportion of the galaxy flux present in the residual images after sub¬ 
tracting the models. In general, cDs have larger Re and RFF values than ellipticals. 
Therefore we find, in a statistically robust way, a boundary separating cD and non-cD 
BCGs in this parameter space. BCGs with cD morphology can be selected with rea¬ 
sonably high completeness (~ 75%) and low contamination 20%). This automatic 
and objective technique can be applied to any current or future BCG sample with 
good quality images. 

Key words: galaxies: clusters: general — galaxies: elliptical and lenticular, cD — 
galaxies: structure 


1 INTRODUCTION 


The brightest cluster galaxies (BCGs) are the most lumi¬ 
nous and massive galaxies in today’s universe. Their stellar 
masses reach beyond ~ IO'^^Mq, and they reside at the bot¬ 
tom of the gravitational potential well of galaxy clusters and 
groups. Their formation and e volution relate clos ely to the 
evolution of the host clusters (IWhilev et abll^OSh and fur- 
ther tie to the hist ory of large-scale structures in universe 
llConrov et al.ll2007^. BCGs are typic ally classified as ellipti¬ 
cal galaxies ( Lauer fc Postmanlll992h , but a fraction of them 
possess an extended, low surface brightness envelope around 


* E-mail: : ppxdzl@nottlnghcmi.ac.uk 
t E-mail: : alfonso.aragon@nottinghain.ac.uk 
J E-mail: : christopher.conseliceanottlnghani.ac.uk 


the central region. These are referr ed to as cD galaxies (e.g. 
lDresslerlll984 lOeeerle fc HillllioOlh . 


The surface brightness profile of elliptical galaxies was 
origi nally modelled using th e empirical ^3 Vaucou leurs 
law (Ide Vaucoule ur^ ligdal l. However, iLuggeil (Il984l ~) and 
ISchombertI ( 1986h showed that the A rnodel cannot prop¬ 
erly describe the flux excess at large radii for most ellipti¬ 
cal galaxies, and an additional para meter n was i ntroduced 
in the so-called Sersic law (ISersid [196311 . For the 

most massive early-type galaxies, however, a single Sersic 
profile still doe s not reproduce their luminosity distribu¬ 
tion accurately. iGonzalez et al.l ((200^ found that a sam¬ 
ple of 30 BCGs were best fitted using a double ^0 
Vauc ouleurs profile rather than a single Sersic law. Further¬ 
more, (DonzeUTeU^ (1201 il l suggested that a two-component 
model with an inner Sersic and an outer exponential profile 
is required to properly decompose the light distribution of 
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~ 48% of the BCGs in the ir 430 galaixy samp le. A similar 
conclusion was obtained bv ISeigar et al.l ll2007l '). 

The light profiles of BCGs need to be explained by any 
successful model of galaxy formation and evolution. In hier¬ 
archical models of str ncture formation, a tw o-phase scenario 
is currently favoured. iHopkins et al.l (l2009l l proposed that a 
early central starburst could give rise to the bulge (ellipti¬ 
cal) component of these galaxies, while the outer envelope 
was subsequently formed by the violent relaxation of stars 
originating in galaxies wh i ch me rged with the central galaxy . 
Alternatively. fOser et al.l ll2010l ') and ijohansson et al.l (l20ld l 
suggested that intense dissipational processes such as cold 
accretion or gas-rich mergers could rapidly build up an ini¬ 
tially compact progenitor and, after the star formation is 
quenched, a second phase of slower, more protracted evolu¬ 
tion is dominated by non-dissipational processes such as dry 
minor mergers to form the low-surface-brightness outskirts. 

To shed light on the mechanism(s) leading to the forma¬ 
tion of BGGs, especially of cD galaxies, we need to answer 
questions such as: are elliptical and cD BGGs two clearly 
distinct and separated classes of galaxies? if so, are elliptical 
and cD BCGs formed by different processes or in different 
environments? are the extended envelopes of cD galaxies in¬ 
trinsically different structures which formed separately from 
the central bulge? To help answer these questions, in this 
paper we explore statistically how the visual classification 
of BCGs into different morphological classes (e.g., elliptical, 
cD; here referred to as “morphology”), relates to the quanti¬ 
tative structural properties of their light profiles (e.g., effec¬ 
tive radius Re, Sersic-index n; genericafly called “structure” 
in this paper). Moreover, finding an automatic and objective 
way to select cD BCGs is nontriv ial for the future databases 
and study. Recent study such as iLiu et aP ([200^ ide ntified 
cD BCGs by Petrosian parameter profiles l|Petrosianlll976li . 
but their method does not give an unambiguous criterion to 
separate cD galaxies from non-cD BCGs. 

In thi s paper, we visually-c l assify 625 BCGs from the 
sample of Ivon der Linden et al.l (l2007l . hereafter L07) and 
ht accurate models to their light profiles. We find clear links 
between the visual morphologies and the structural param¬ 
eters of BCGs, and these allow us to develop a quantitative 
and objective method to separate cDs galaxies from ellipti¬ 
cals BCGs. In a later paper (Zhao et ah, in preparation) 
we will study how the visual morphology and structural 
properties of BCGs correlate with their intrinsic properties 
(stellar masses) and their environment (cluster mass and 
galaxy density), and explore the implications that such cor¬ 
relations have for the formation mechanisms and histories 
of cDs/BCGs. 

The paper is organized as follows. In lj2] we introduce 
the BCG samples and the visual morphological classification 
of the BCGs. In Sj3]we describe the light distribution models 
and the fitting methods we use, and discuss how the results 
are affected by sky-subtraction uncertainties. This section 
also presents a quantitative evaluation of the quality of the 
fits. In m we present the structural properties of the BCGs 
in the sample. In an we introduce an objective diagnostic 
to separate cDs from non-cD BCGs using quantitative in¬ 
formation from their light profiles. We summarise our main 
conclusions in aH 


2 DATA 

2.1 BCG Catalogue and Images 

To study the structural properties of BCGs in galaxy groups 
and clusters, we use the BCG catalogue published by L07. 
The groups and clusters that host these B CGs come from 
the C4 cluster catalogue ijMiller e t a .Il20 05fl extracted from 
the Sloan Digital Sky Survey (SDSS; York et ^l2000l 'l third 
data release spectroscopic sample. The cluster-finding al¬ 
gorithm used to build the C4 catalogue identifies clusters 
as over-densities in a seven-dimensional parameter space 
of position, redshift and colour, minimising projection ef¬ 
fects. The C4 catalogue gives a very clean widely-used clus¬ 
ter sample which is well supported by simulations. Based 
on the C4 catalogue, L07 restricted their cluster sample 
to the 0.02 ^ 2 < 0.10 redshift range to avoid prob¬ 
lems related to the 55 arcsec “fiber collision” region of 
SDSS. Within each cluster, L07 applied an improved algo¬ 
rithm to identify the BCG as the galaxy being closest to 
the deepest point of the p otential well of the cluster (see 
Ivon der Linden et al.ll20(^ for a detailed discussion of this 
identification), and developed an iterative algorithm to mea¬ 
sure the cluster velocity dispersion ar 200 within the virial ra¬ 
dius i?20(ll|. The catalogue created by L07 contains 625 BCGs 
in galaxy groups and clusters with redshifts 0.02 ^ 2 ^ 0.10 
spanning a wide range of cluster velocity dispersions, from 
galaxy groups (o'r 200 ^ 200 km/s) to very massive clusters 
(crr 200 ~ 1000 km/s). 75% of the BGGs in L07 are in dark 
matter halos with CTr 200 ^ 309 km/s, where the completeness 
of the halos identihed by the G4 algorithm is expected to be 
above 50%. Obviously, for larger halo masses the complete¬ 
ness is higher. 

The images we use to classify the BCGs and analyse 
their structural properties come from the SDSS Seventh 
Data Release (DR7) r-band images. We also use SDSS-DR7 
p-band images of the BGGs in Section 14.11 The BCG cata¬ 
logue used in this paper together with their main properties 
are presented in appendix A. 

2.2 Visual Classification 

The 625 BCGs in L07 sample were visually classihed by 
careful inspection of the SDSS images. BCGs were displayed 
using a logarithmic scale between the sky level and the peak 
of the surface brightness distribution. The contrast was ad¬ 
justed manually to ensure that the low-surface-brightness 
envelopes were revealed if present. cD galaxies are identified 
by a visually extended envelope, while the envelope is not 
visible in our elliptical BGGs. Finally the BGGs were classi¬ 
hed into three main types: 414 cDs, including pure cD (356), 
cD/E (53) and cD/SO (5); 155 ellipticals, including pure E 
(80), E/cD (72), and E/SO (3); 46 disk galaxies, containing 
spirals (24) and SOs (22). The main morphological classes of 
BGGs are illustrated in Figure [T] There are also 10 BGGs 
undergoing major mergers, but we will not discuss them in 
this paper in any detail. 

Over half of the BCGs in the sample are classihed as 
cDs. Separating cD BGGs and non-cD elliptical BCGs is a 

^ R 200 is the radius within which the average mass density is 
200pc, where pc is the critical density of the universe. 
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Figure 1. Examples of the main morphological classes of BCGs in our sample (cD, cD/E, E/cD, E, cD/SO, E/SO, SO, Spiral) illustrating 
the gradual transition between classes. The images are displayed using a logarithmic surface-brightness scale. 


very hard problem since there is no sharp distinction be¬ 
tween these two classes (e.g., IPatel et al.[ l200^ : iLiu et ahl 
l2008l ). Detecting the extended stellar envelope that char¬ 
acterises cD galaxies depends not only on its dominance, 
but also on the quality and depth of the images, and on 
the details of the method(s) employed. We used intermedi¬ 
ate classes such as cD/E (probably a cD, but could be E) 
and E/cD (probably E, but could be cD) to account for the 
uncertainty inherent in the visual classification. 

Our careful inspection of the images clearly reveals that 
there is a wide range in the brightness and extent of the 
envelopes. There seems to be a continuous distribution in 
the envelope properties, ranging from undetected (pure E 
class) to clearly detected (pure cD class), with the interme¬ 
diate classes (E/cD and cD/E) showing increasing degrees 
of envelope presence. This continuous distribution in enve¬ 
lope detectability will also be made evident in the structural 
analysis carried out later in this paper. The classification 
we present here does not intend to be a definitive one since 
such a thing is probably unachievable. Our aim is to ob¬ 
tain a homogeneous and systematic visual classification of 
the BCGs and then study how such classification correlates 
with quantitative and objective structural properties of the 
BCGs. The visual morphological types of all the galaxies in 
the sample are presented in Appendix A. 

We checked the effect that the redshift of BGGs may 
have on the visual classification. cDs might be mistakenly 
identified as elliptical if they are more distant since the ex¬ 
tended low-surface-brightness envelope may be harder to re¬ 
solve at higher redshifts. Figure [2] illustrates the redshift 
distribution of the three main types. cD galaxies generally 
share the same redshift distribution with elliptical BGGs, 
especially at z ^ 0.05. At g < 0.05 we identify a slightly 
higher proportion (by ~ 10%) of cD galaxies. However, if 
we compare the structural properties of cD and elliptical 
BCGs which are at z ^ 0.05, the results we obtain do not 
significantly differ from those using the full-redshift sample. 
As an additional check, we artificially redshifted some of the 



Figure 2. Redshift distribution for BCGs with different morpho¬ 
logical types. The red solid line corresponds to cD BCGs, the 
green dashed line to ellipticals, and the blue dotted line to disk 
(spiral and SO) BCGs. A Kolmogorov-Smirnov test indicates that 
the redshift distributions of cD and elliptical BCGs are only dif¬ 
ferent at the ~ 2.4 (t level. cD galaxies share the same redshift 
distribution with elliptical BCGs at z ^ 0.05, while there are 
proportionally ~ 10% more cD galaxies at z < 0.05. 


lowest redshift galaxies (z ~ 0.02-0.03) to z = 0.1, the high¬ 
est redshift of the sample, taking into account cosmological 
effects such as surface-brightness dimming. Because the red¬ 
shift range of the BCGs we study is very narrow, the effect 
on the images is minimal and does not have any significant 
impact on the visual classification. We are therefore confi¬ 
dent that our visual classification is robust and that in the 
relatively narrow redshift range explored here any putative 
redshift-related biases will not affect our results. 
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3 QUANTITATIVE CHARACTERISATION OF 

BCG STRUCTURE 

The surface brightness profiles of galaxies provide valuable 
information on their structure and clues to their formation. 
It has become customary to fit the radial surface bright¬ 
ness distribution using theoretical functions which have pa¬ 
rameters that include a measurement of size (e.g., half-light 
radius or scale length), a characteristic surface brightness, 
and other parameter(s) describing the shape and proper¬ 
ties of the surface brightne ss profiles. In this paper we use 
GALFIT (IPeng et al.ll200 to fit the 2-D luminosity profile 
of each BCG using two parametric models, and thus deter¬ 
mine the best-fitted parameters of each model. GALFIT can 
simultaneously fit model profiles to several galaxies in one 
image, which is particularly important for BCGs since they 
usually inhabit very dense environments. In this way, the 
light contamination from nearby galaxies can be accounted 
for appropriately. 

We explore two models to represent the luminosity pro¬ 
file of the BCGs. A model commonly used to fit a vari- 
ety of galaxy li g ht pr ofiles is the general ization of th e 
Ide VaucouleursI (|l948h law introduced bv ISersij (llQb.lh . The 
Sersic model has the form 

7(r) = 7eexp{-6[(r/re)^'^" - 1]}, (1) 

where 7(r) is the intensity at distance r from the centre, 
Re, the effective radius, is the radius that encloses half 
of the total luminosity, 7e is the intensity at Re, n is the 
Sersic index repre senting concentration, and fo ~ 2n — 0.33 
dCaon et al.l ll993l L The Sersic function provides a good 
model for galaxy bulges and massive elliptical galaxies. Since 
BCGs are mostly early-type galaxies, it is reasonable to fit 
their structure with single Sersic models first. Subsequently, 
in order to explore the complexity introduced by the ex¬ 
tended envelopes of cD galaxies, we will also fit the light 
profile of BCGs adding an additional exponential compo¬ 
nent to the Sersic profile. Adding this exponential compo¬ 
nent is the simplest way to describe the “extra-light” from 
the extended envelope. Note that the exponential profile 
7(r) = 7o exp(—r/rs) is just a Sersic model with n = 1. The 
models assume that the isophotes have elliptical shapes, and 
the ellipticity and orientation of each model component are 
parameters determined in the fitting process. 

In order to run GALFIT, we require a postage stamp 
image for each BCG with appropriate size to measure its 
structure over the full extent of the object, a mask image 
with the same size as the stamp image, an initial guess for 
the fitting parameters, an estimate of the background sky 
level, and a point spread function (PSF). Details on how 
these ingredients are produced and the fitting procedures 
are given below. 


3.1 Pipeline for One-Component Fits: Modified 
GALAPAGOS 

We run GALFIT using the GALAPAGOS pipeline 
llBarden et al.|[^012h . GALAPAGOS has been successfully 
applied to a wide variety of ground- and space-based image s 
dHaufiler et al.l 20131 : IVika et ^ l2013l : IVulcani et al.l l2014h . 
IGuo et ^ ( 2 OO 9 II specifically tested a modified version of 
GALAPAGOS to fit the central galaxies of local clusters 


using SDSS r-band images. We applied the same modified 
version of GALAPAGOS to fit the SDSS r-band images of 
the BCGs in our sample. The starting point are SDSS im¬ 
ages with a size of 2047 x 1488 pixels. For each BGG, the 
pipeline carries out four main tasks before running GALFIT 
itself: (i) detection of all the sources present in the image; 
(ii) cutting out the appropriate postage stamp and prepar¬ 
ing the mask image; (iii) estimation of the sky background; 
(iv) preparation of the input file for GALFIT. After com¬ 
pleting these tasks, GALAPAGOS will run GALFIT using 
the appropriate images and input parameters. We describe 

now these tasks in detail. _ 

_ (i) Source Detection: SExtractor llBertin fc Arnorit3 

Il996h is used to detect galaxies in the SDSS images. A 
set of configuration parameters defines how SExtractor de¬ 
tects source s. The values of the SExtractor input parame¬ 
ters follows IGuo et al.l (l2009lL DETEGT_MINAREA = 25, 
DETEGT.THRESH = 3.0, and DEBLEND_MINCONT = 
0.003. This set of parameters were tested to perform well 
on SDSS r-band images so that the bright and extended 
BCGs were isolated from other sources without artificially 
deblending them into multiple components. SExtractor also 
provides estimates of several properties for the target BCGs 
and nearby objects such as their magnitude, size, axis ratio 
and position angle. These values are used to calculate the 
initial guesses of the model parameters that are needed as 
inputs by GALFIT. 

(ii) Postage Stamp creation: GALAPAGOS cuts out a 
rectangular postage stamp centred on the target BCG which 
will be used by GALFIT as input image. We define the 
“Kron ellipse” for a galaxy im^e as an ellipse whose semi¬ 
major axis is the Kron radiuul (7?kron), with the elliptic¬ 
ity and orientation determined by SExtractor. The postage 
stamp size is determined in such a way that it will fully con¬ 
tain an ellipse 3.5 times larger than the Kron ellipse, i.e., 
its semi-major axis is 3.57?kron, and has the same ellipticity 
and orientation. The 3.5 factor represents a compromise be¬ 
tween computational speed and ensuring that virtually all 
the BGG’s light is included in the postage stamp. At this 
stage, a mask image is also created, identifying and mask¬ 
ing out all pixels belonging to objects in the postage stamp 
which will not be simultaneously fitted by GALFIT. The 
aim is to reduce the computational time by excluding ob¬ 
jects too far from the BCG or to o faint to have any sig¬ 
nificant effect on the fit. Following iBarden et al.l (l2012h . an 
“exclusion ellipse” is defined for each galaxy with a semi¬ 
major axis 1.57?kron +20 pixels, and the same ellipticity and 
orientation as the Kron ellipse. GALAPAGOS masks out all 
objects whose exclusion ellipse does not overlap with the ex¬ 
clusion ellipse of the target BGG. These objects are deemed 
to be too far away from the BGG to require simultaneous 
fitting. Furthermore, all objects more than 2.5 magnitudes 
fainter than the BCG are also masked out since they are too 
faint to affect the BCG fit. The pixels that belong to these 
objects according to the SExtractor segmentation maps are 
masked out and excluded from the fits. All the remaining ob- 


^ In this paper we use the following definition of “Kron radius”: 
7?krot i = 2.5ri, where r\ is the first mome nt of the light distribu¬ 
tion dKro r] |l980l : [Bertin Arnoutg||l996h . For an elliptical light 
distribution, this is, strictly speaking, the semi-major axis. 
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jects will be simultaneously fitted by GALFIT at the same 
time as the BCG. For a detailed description of this process 
and a justification of the parameter choice see lBarden et al.l 

ll2012tl . 

(in) Sky Estimation: Accurate estimates of the sky 
background level is crucial when fitting galaxy profiles, par¬ 
ticularly when interested in the low-surface-brightness outer 
regions. Overestimating the sky level will result in the under¬ 
estimation of the galaxy flux, size, and Sersic index n, and 
vice-versa. GALAPAGOS uses a flux growth curve method 
to robustly estimate the local sky background around the 
target galaxy. SDSS DR7 also provides a global sky value 
for the whole 2047 x 1488 image frame and local sky values 
for each galaxy. The SDSS PHOTO pipeline estimates the 
sky background using the median flux of all the pixels in 
the image after 2.33cr-clipping. However, according to the 
SDSS-III website, the version of PHOTO used in DR7 and 
earlier data releases tended to overestimate both the global 
and local sky values. The sky measurement is improved by 
SDSS-III in later data releases, but since we use the images 
from DR7 we cannot use t he SD SS sky value with enough 
confidence. iHaussler et aP |2003) demonstrated that the sky 
measurement that GALAPAGOS produces is highly reliable 
for single-band fits because it takes into account the effect 
of all the objects in the image. Therefore, in this study we 
use the local sky background estimated by GALAPAGOS. 
The accurate sky measurement provided by GALAPAGOS 
indicates that we can reach a surface brightness limit in the 
r-band of ~ 27 mag/arcsec^. This is deep enough to study 
the faint extended structures of BCGs. For each galaxy, its 
local sky background is included in the GALFIT input file 
and is fixed during the fitting procedure. Given the impor¬ 
tance of accurate sky subtraction, in Section 13.21 we will 
carry out an explicit comparison of our results using SDSS 
and GALAPAGOS sky estimates. 

(iv) GALFIT Input: GALAPAGOS produces an input 
file which includes initial guesses for the fitting parameters 
based on the SExtractor output. As mentioned above, all 
objects which are not masked out are fitted simultaneously 
using a Sersic model. The initial-guess model parameters for 
these nearby companions are also determined from SExtrac¬ 
tor. In order to obtain reasonable results, we impose some 
constraints on the acceptable model parameter range. Our 
constrains o n position, magnitude , axis ratio and position 
angle follow IHaussler et al.l (|2007l i. Additionally, the half- 
light radius Re is constrained within 0.3 ^ i?e 800 pixels. 
This prevents the code from yielding unreasonably small 
or large sizes. Since the pixel size of the SDSS images is 
0.396 arcsec, Re is constrained to be larger than 0.12arcsec, 
which is much smaller than the PSF, and smaller than half 
the size of the original input images, reasonable for the range 
of redshifts explored. In the original GALAPAGOS pipeline, 
the constraint on the Sersic index is 0.2 ^ n ^ 8. These 
are reasonably conservative limits, since normal galaxies 
with n > 8 are rarely seen and are often associated with 
poor model fits. However, some studies have shown that 
very luminous ellipt ical galaxies with n > 8 do exist (e.g., 
iGraham et al.ll2005l i. therefore for the target BCGs we al¬ 
low n to be as large as 14 to keep the fits as free as possible. 
For the companion galaxies, which are fitted simultaneously, 
we still keep the constraint 0.2 n 8. The final ingre¬ 
dient needed by GALFIT is a PSF image appropriate for 



sky,dss-skyg,,i,, (ADU) 

Figure 3. Distribution of the difference between the SDSS DR7 
global sky and the GALAPAGOS-measured sky values. In gen¬ 
eral, SDSS overestimates the sky background. The average sky 
value measured by GALAPAGOS in the SDSS r-band BCG im¬ 
ages is 140.8 ADU, corresponding to a surface brightness of ~ 20.9 
mag/arcsec^. 

each BCG. These are extracted from the SDSS DR7 data 
product^ according to the photometric band used and the 
position of the BCG on the SDSS image. 

3.2 Effect of the Sky Background Subtraction: 

Comparing SDSS and GALAPAGOS Sky 

Estimates 

As described in Section [Ml in this study we rely on the sky 
measurements provided by GALAPAGOS. However, it is im¬ 
portant to test the effect that the choice of sky background 
has on our results. We do this by comparing the htted Sersic 
model parameters n and Re using the GALAPAGOS and 
SDSS sky estimates. As mentioned before, SDSS DR7 pro¬ 
vides a global sky value for the w hole 2047 x 1488 image 
and local sky values for each galaxy. I Guo et ^ ll2009l 'l found 
that the local background estimates are generally larger than 
the global ones due to contamination from the outskirts of 
extended and bright sources, making them unreliable. We 
therefore restrict our comparison to the global SDSS sky 
values. We ht the BCG light profiles twice using exactly the 
same procedure and input parameters (see 1)3.111 but chang¬ 
ing only the sky background estimates. The hrst set of fits 
use the GALAPAGOS-determined sky values, while the sec¬ 
ond set use the SDSS DR7 global ones. 

Figure [3] shows the distribution of the difference be¬ 
tween the SDSS DR7 global sky and the sky measured by 
GALAPAGOS. It is clear that the SDSS global sky is gen¬ 
erally larger than the local sky from GALAPAGOS. The 
effect from different sky values on the best-htted structural 
parameters (Sersic index n and effective radius Re) is shown 

® http://www.sdss.org/DR7/products/images/read_psf.html 
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Figure 4. Comparison on the best-fit n and Re from single Sersic 
models using the SDSS and GALAPAGOS-measured sky esti¬ 
mates. Solid and open red circles correspond to pure cD and other 
cD galaxies (cD/E and cD/SO) respectively; solid and open green 
diamonds correspond to pure and other (E/cD and E/SO) ellip¬ 
tical BGGs respectively; solid blue triangles represent SOs and 
open ones are spirals. It shows that the SDSS overestimation of 
the global sky result in the values of rigUgg and rg^sdss being smaller 
than the corresponding GALAPAGOS ones. Moreover, the effect 
is more serious for the BGGs with large n and Re which are mostly 
cDs. 


in Figure [H It is clear that the SDSS larger sky values re¬ 
sult in the values of nsdss and re.sdss being smaller than 
the corresponding GALAPAGOS ones. The effect becomes 
more severe for those BGGs with large n and Re, most of 
which are cD galaxies. This means the overestimated sky 
values would particularly affect the measurements on the 
low-surface-brightness envelopes of cD galaxies. Although 
it is difficult to know a priori which the true value of the 
sky background is, based on the fact that the SDSS-III pro- 
vides evidence t hat D R7 sky values are overestimated while 
lllaussler et al.l ll2007l ) showed reasonable proof of the re¬ 
liability of the GALAPAGOS sky measurements, in what 
follows we will therefore trust and use the GALAPAGOS- 
determined sky values. 


3.3 Two-Component Fits 

Although the light profiles of many early-type galaxies can 
be reproduces reasonably well with single Sersic models, the 
extended envelopes of cD galaxies may require an additional 
component. We therefore fitted all the BGGs using a two- 
component model consisting of a Sersic profile plus an ex¬ 
ponential. The input postage stamp, mask image, PSF, and 
sky values required by GALFIT remain the same as for the 
single-Sersic fits. To ensure that we are fitting exactly the 
same light distribution, the location of the centre of the BGG 
is fixed to the X and Y coordinates determined in the single 


fit, and we also force the initial guesses of the model pa¬ 
rameters to be the single-component fit results. The BCG 
companions are simultaneously fitted still with single-Sersic 
profiles but with initial-guess parameters determined by the 
single profile fits. 


3.4 Residual Flux Fraction and Reduced 

Although the models we are fitting are generally reasonably 
good descriptions of the BCG light profiles, real galaxies 
can be more complicated, with additional features and struc¬ 
tures such as star-forming regions, spiral arms, and extended 
halos. It is therefore desirable to quantify how good the fits 
are and what residuals remain after subtracting the best-fit 
models. A visual inspection of the residual images can gen¬ 
erally give a good feel for how good a fit is, and sometimes 
tell us whether an additional component or components are 
required. However, more quantitative, repeatable and objec¬ 
tive di agnostics are also needed. The residual flux fraction 
1RFF: IHovos et ^120111 ) provides one such diagnostic. It is 
defined as 


RFF = 


E. 


iJeA 


\Ii,j - - 0-8 X 


bkg 




( 2 ) 


where A is the particular aperture used to calculate RFF. 
Within A, lij is the original flux of pixel (i.j), is 

the model flux created by GALFIT, and is the rms of 
the background. RFF measures the fraction of the signal 
contained in the residual image that cannot be explained 
by background noise. The 0.8 factor ensures that the expec¬ 
tation value of the RFF for a purely Gaussian nois e error 
image of constant variance is 0.0. See IHovos et al.l (1201 ll) 
for details. Obviously, this diagnostic can be applied to 
both single-Sersic and two-component profiles, or any other 
model. The aperture A we use to calculate RFF is the “Kron 
ellipse” defined in Section o (an ellipse with semi-major 
axis Rkron and the ellipticity and orientation determined 
by SExtractor for the BCG). the denominator of 

Equation is computed as the total BCG flux contained 
inside the Kron ellipse, which is one of the SExtractor out¬ 
puts, and therefore independent of the model fit. 

Since BGGs usually reside in dense environments, some¬ 
times there are some faint nearby objects contained within 
the Kron ellipse that have not been fitted by GALFIT (those 
more than 2.5 mag fainter than the BCG, see >13.111 . These 
objects will be present in the residual image. Moreover, 
brighter companions that have been simultaneously fitted 
may also leave some residuals due to inaccuracies in their 
fits. Therefore, even if the BCG light distribution has been 
accurately fitted, RFF can be affected by the residuals from 
the companion galaxies, failing to provide an accurate mea¬ 
sure of the quality of the fit. To minimise the effect from 
companion galaxies on RFF, we mask out the pixels belong¬ 
ing to all companions within the Kron ellipse using SExtrac¬ 
tor segmentation maps. The RFF will therefore measure the 
residuals from the BCG fit alone, excluding, as far as possi¬ 
ble, those belonging to nearby galaxies. 

An additional measurement of the fit accuracy is the 
reduced x^i which is minimised by GALFIT when finding 
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the best-fit models. It is defined as 


X 


2 

u 


1 

Ndof 


E 


{L.j 


rmodel\2 

-^i,j) 



(3) 


where A is the aperture used to calculate -^dof is the 
number of degrees of freedom in the fit, lij is the original 
image flux of pixel represents, for each pixel, the 

sum of the flux of the models fitted to all the galaxies in the 
aperture, and Uij is the noise corresponding to pixel {i,j). 
This noise is calculated by GALFIT taking into account the 
contributio n of the Poisson errors and the read-out-noise of 
the image dPeng et al.ll200j i. 

Similarly to RFF, Xu also measures the deviation of the 
fitted model from the original light distribution. The value 
of Xv that GALFIT minimises to Hnd the best-fit model 
is calculated over the whole postage stamp, and includes 
contributions from ah the objects htted. To make sure that 
we only take into account the contribution to Xv from the 
BCG fit, we calculate it within the Kron ellipse of the BCG, 
masking out the nearby objects as we did when calculating 
RFF. 

The choice of aperture (Kron ellipse with semi-major 
axis of i?kron) Over which we evaluate RFF and xf represents 
a good compromise between covering a large fraction of the 
galaxy light while minimising the impact of close compan¬ 
ions. We carried out several tests to evaluate the sensitivity 
of our results to the changes in aperture size. If we reduce 
the semimajor axis of the aperture by 20% or more we lose 
signihcant information on the extended halo of BCGs, which 
we must avoid. If we increase the semimajor axis of the aper¬ 
ture by 20% or more, we potentially increase the sensitivity 
to the galaxy halos but in the crowded central cluster regions 
contamination from companion galaxies becomes a serious 
problem, generally increasing RFF and xf. Changes in the 
aperture semimajor axis within ±20% would have no effect 
on the conclusions of this paper. 


3.5 Evaluating One-Component and 
Two-Component Fits 


residuals often show excess light at large radii which can be 
identified as an exponential component or halo. Addition¬ 
ally, the fit to these galaxies visibly improves when using 
a two-component model. With these criteria we confidently 
identify 53 1C BCGs and 25 2C BCGs. Since we want to 
test the sensitivity of RFF and x?i we concentrate for now 
on this small but robust subsample. The rest of the BCGs 
(537) cannot be confidently classified into 1C or 2C BCGs 
because it is too hard to tell visually due to the residuals 
containing significant structures which cannot be accurately 
fitted by such simple models. 

Figure [5] presents a comparison of the RFF and xL val¬ 
ues for the one- and two-component fits of the 53 1C BCGs 
and 25 2C BCGs. For 1C BCGs, the RFF and xS distribu¬ 
tions of one- and two-component fits are virtually indistin¬ 
guishable. Neither RFF nor xt improve significantly when 
the second component is added. However, RFF and x2 are 
significantly smaller for the two-component fits of 2C BCGs. 
It is clear therefore that the quantitative information that 
RFF and xf provide agrees very well with the visual as¬ 
sessments of the fits. Both RFF and xf are sensitive to 
changes in the residuals, but RFF appears to be more sen¬ 
sitive. As shown in the bottom panels of Figure |6l the im¬ 
provement in the two-component fit for 2C BCGs is around 
40%-60% when measured by RFF, while it is only ~ 20% 
when measured by xS- A further useful piece of information 
obtained from this test is that the typical values of \ogRFF 
and logxJ for fits deemed to be good by visual inspection are 
log7?FF~ -l-Tio.Se. and of logxl'^ 0.0421°;“^ (median 
+/— 1st and 3rd quartiles of the parameter distributions). 

As mentioned before, the majority of the BCGs cannot 
be visually classified into 1C or 2C BCGs with high certainty 
because their light distributions are too complex to be ac¬ 
curately represented by such simple models. Nevertheless, 
we can use the quantitative information provided by RFF 
and Xv to gauge to what extent the BCGs are better fit by 
a two-component model than by a one-component model. 
This will be discussed later. 


Since RFF and xS can quantify the residual images after 
subtracting the model fits, we attempt to use them to assess 
whether a one-component (Sersic) fit or a two-component 
(Sersic±Exponential) fit is more appropriate to describe the 
light profile of individual BCGs. In order to do this, we first 
evaluate the effectiveness of RFF and xt s-t quantifying the 
goodness-of-fit. We visually examine the fits and residuals 
obtained from both one- and two-component models for all 
the BCGs in our sample. In some cases, two of which are 
illustrated in Figure [5l it is obvious which model is clearly 
favoured. 

For those BCGs where such a clear distinction can con¬ 
fidently be made, we classify them into what we call 1C (one- 
component) BCGs and 2C (two-component) BCGs. Explic¬ 
itly, 1C BCGs (e.g., galaxy 1 in the top panel of Figure (Sj 
are those for which a one-component Sersic model represents 
their light distribution very well, and therefore the residuals 
left are small and show no significant visible structure. For 
these galaxies, adding a second component does not visibly 
improve the residuals. Conversely, 2C BCGs (e.g., galaxy 2 
in the bottom panel of Figure[5]) are not well fitted by a one- 
component model, and the residuals are significant. These 


We would like to point out that this is the first time that 
the residual flux is calculated considering only the contribu¬ 
tion of the target galaxies when estimating both RFF and 
xt, explicitly excluding the contribution due to the com¬ 
panion galaxies. For instance, IHovos et al.l (1201 il l also used 
RFF to evaluate the goodness-of-fit, but they measured the 
residuals over all pixels within a specific area around the 
target galaxies, without excluding nearby companions. Sim¬ 
ilarly, the xf values from GALFIT have also been applied 
to ev aluate which fitting model is better (e.g., iBruce et al.l 
I2OI2II . but the effect of nearby objects on the xf values was 
also overlooked. Using the 2C BCG sample, we assessed the 
importance of this improvement. If the RFF and xS are cal¬ 
culated considering the residuals in all the pixels inside the 
relevant aperture, the RFF and xJ distributions for the two- 
component fits of 2C BCGs cannot be distinguished from 
the one-component results. The effect of the contribution 
to the residuals from companion galaxies is so severe that 
it renders such a comparison useless. Our method therefore 
represents a significant step forward. It is extremely impor¬ 
tant to exclude the contibution of the companion galaxies 
when calculating RFF and xf in this kind of analysis. 
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Figure 5. Example of one-component (Sersic) fits and two-component (Sersic-I-Exponential) fits for 1C and 2C BCGs, respectively. 
From left to right, the panels show the original image, the one-component model, the residuals after subtracting the one-component 
fit, the two-component model, and the residuals after subtracting the two-component fit. The upper panels show a 1C BCG where a 
one-component fit does a good job and adding a second component does not visibly improve the residuals. The lower panels show a 
2G BGG, where the one-component residual exhibits clear excess light at large radii, suggesting that a second component is necessary. 
Indeed, the two-component residual is much better for this BGG. 


4 STRUCTURAL PROPERTIES OF BCGS 

Our morphologically-classified BCGs provide a large sam¬ 
ple to statistically study their structural properties and link 
them to their morphological properties. In what follows we 
consider the three main morphological classes of BCGs: cDs 
(including all BCGs classified as pure cD, cD/E and cD/SO); 
ellipticals (including pure E, E/cD and E/SO) and disk (spi¬ 
ral and SO) BCGs. The 10 BCGs classified as mergers are 
excluded (see Section [2. 2 1 for details). We decided to include 
the galaxies with “uncertain” morphologies (such as cD/E 
and E/cD) in our analysis to reflect the difficulties involved 
in visual classification. However, to ensure the robustness of 
our analysis, at every stage we have checked that consider¬ 
ing only “pure” cD and elliptical BCGs (i.e., excluding the 
cD/E, cD/SO, E/cD and E/SO classes) would not change our 
conclusions. 

Since most BCGs are early-type galaxies, we will 
first consider and discuss single Sersic models when fit¬ 
ting their SDSS r-band images. We will subsequently use 
Sersic-I-Exponential models to see whether the fits are im¬ 
proved. But before embarking in the analysis of the param¬ 
eters derived from these model fits, we first evaluate their 
uncertainties. 


4.1 Structural Parameter Uncertainties 

The parameter uncertainties that GALFIT reports are cal¬ 
culated using the covariance matrix derived from the Hes¬ 
sian matrix computed by th e Levenberg-Marq uardt algo¬ 
rithm that the program uses (IPeng et al.N201oh . These for¬ 
mal uncertainties are only meaningful when the model pro¬ 
vides a good fit to the image, in which case the fluctuations 
in the residual image are only due to Poisson noise. However, 
for real galaxy images the residual images contain not only 
Poissonian noise, but also systematics from non-stochastic 
and stochastic factors due to additional components not in¬ 
cluded in the fitting function (e.g., spiral arms, star-forming 
regions), asymmetries, shape mismatch, flat-fielding errors 


and so on. These non-random factors usually dominate the 
uncertainty of the parameters, and the uncertainties inferred 
from the covarian ce matrices are only lower-limit estimates 
dPeng et al.ll201ol l. Therefore, if we rely on the errors re¬ 
ported by GALFIT the uncertainties in the structural pa¬ 
rameters of the BCGs could be severely underestimated. In¬ 
deed, these formal errors seem unrealistically small: typical 
GALFIT uncertainties for Re and n are only ~ 1-2%. A 
more robust and realistic way of determining these uncer¬ 
tainties is clearly needed. 

We have measured the structural parameters of the 
BCGs in our sample using the SDSS r-band images. Inde¬ 
pendent measurements can also be obtained using the SDSS 
p-band images. In principle, the structural parameters could 
be wavelength-dependent. However, the g — r colours of mas¬ 
sive early-type galaxies with old stellar populations are quite 
spatially un iform and do not cha nge much from galaxy-to- 
galaxy (e.g.. lFukugita et al.lll99^ L Furthermore, morpholog¬ 
ical fc-corrections are neglig ible for early-type galaxi es be¬ 
tween these two bands fe.g.. lTavlor-Mager et al.]|2007^ . so it 
is reasonable to expect that the intrinsic structural param¬ 
eters will not change much between g and r band. There¬ 
fore, any differences in the measured parameters between 
these two bands should be largely dominated by measure¬ 
ment errors. Moreover, if there are significant wavelength- 
dependent differences in the measured parameters that are 
driven by real physical differences, it is reasonable to ex¬ 
pect that these may correlate with other galaxy properties 
such as their colour, morphology, redshift, cluster velocity 
dispersion, etc. No such correlations were found, so we are 
confident that the intrinsic differences are not significant in 
these two bands. 

We use GALAPAGOS to fit the SDSS p-band images 
of the BGGs in our sample in exactly the same way as we 
did for the r-band images. Figure [3 shows a comparison of 
the Re and RFFic values obtained in both bands. Similar 
comparisons were carried out for the rest of the structural 
parameters. The scatter around the 1-to-l relations is due, 
in principle, to both intrinsic wavelength-dependent differ- 
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Figure 7. Comparison of the Re (left panel) and RFFie (middle panel) values obtained in both the SDSS g-band and r-bands. The 
solid lines correspond to the 1-to-l relations. The right panel shows log(i?e,g/fie,r) vs. log{RFFic,g/RFFic,i). The error bars in the 
bottom-right corner are derived from the rms scatter of each parameter. 


ences and measurement errors. Since, as we have argued, 
the intrinsic differences are not expected to be significant 
between these two bands, the measurement errors should 
dominate the scatter. We can thus use this scatter as an 
estimate of realistic, albeit perhaps marginally pessimistic, 
parameter uncertainties. The average errors are S{n) ~ 0.9, 
5(logre) ~ 0.16, and S{logRFFic) ~ 0.13. 

The right-hand panel of Figure [7] shows that the errors 
in Re and RFF are not correlated. This is an important 
point since these two are the main parameters that we will 
use as diagnostics in our analysis in Section [5] 

4.2 Single Sersic Models 

We analyse now the behaviour of four parameters derived 
from the best-fitting single-Sersic models along with the 
morphological classifications. Two of them, the Sersic in¬ 
dex n and the effective radius Fie, provide information on 
the intrinsic properties of the BCGs. The other two, RFF 
and X?, show how well the models fit the real light distribu¬ 
tion of the BCGs and also provide information about their 
detailed structure. The values of these parameters are listed 
in Appendix A. Figure |8] shows the distribution of these pa¬ 
rameters for the three main BCG morphologies. The a value 
in each panel indicates the significance (confidence level) of 
the observed differences between the cD and elliptical BCG 
parameter distributions. These are derived from two-sample 
Kolmogorov-Smirnov tests. 

4-2.1 Sersic Index n 

The Sersic index n measures the concentration of the light 
profile, with larger n corresponding to higher concentration. 
The upper left panel of Figure |S] presents the n distributions 
for the three main BGG morphologies. It is clear that disk 
(spiral and SO) BGGs tend to have smaller values of n, as ex¬ 
pected. However, the n distribution for disk BCGs is skewed 
towards larger values (n > 3) than t hose of the norm al disk 
galaxy population (e.g., n = 2.5 in IShen et akllioo^ ). This 
is because most disk BCGs are early-type bulge-dominated 
spirals and SOs. Elliptical and cD BGGs tend to have larger 
n values (n 4). The n distributions of cD and elliptical 


BCGs are quite similar. A K-S test indicates that the dis¬ 
tributions are not significantly different: the significance of 
any possible difference is just 2.04 ct. 

4-2.2 Effective Radius Re 

The effective radius Re is a measurement of the extent (or 
size) of the light distribution. The upper right panel of Fig¬ 
ure |5] shows the distributions of log Re- Disk BCGs tend to 
have relatively small sizes, and the vast majority of them 
(~ 85%) have Re smaller than ~ 15h“^kpc. About 75% 
of the elliptical BGGs also have Re < 15h“^kpc, while cD 
galaxies tend to be significantly larger. More than 60% of 
cDs have Re > 15h“^kpc. A K-S test demonstrates that 
the difference in Re distributions between cD and elliptical 
BCGs is very significant. This suggests that Re could be a 
good discriminator to separate cD and elliptical BCGs. 

4-2.3 Residual Flux Fraction and Reduced 

The lower left panel of Figure [8] presents the RFFic distri¬ 
butions in a logj^g scale, where RFFic denotes RFF for one- 
component models. The RFFie of disk BGGs has a much 
broader distribution and reaches significantly larger values 
than those of cDs and ellipticals. This reflects the fact that 
a single-Sersic model is not a good representation of the 
light distribution of galaxies with clear disks, spiral arms and 
star-forming regions. Early-type BCGs have smoother light 
distributions that can be reasonably well reproduced with 
a Sersic profile, and their RFFic tend to be smaller. How¬ 
ever, there are statistically significant differences between 
the RFFie distributions of cD and elliptical BCGs. About 
60% of elliptical BGGs have RFFic values in the range cor¬ 
responding to good fits (see Section [3. 5 1 and Eigure[S]l, while 
just ~ 25% of cD galaxies do. This suggests that most ellip¬ 
tical BCGs can be well represented by single Sersic models, 
while most cD galaxies are harder to model with such a sim¬ 
ple profile. Since an extended envelope is a general property 
of cD galaxies, their deviation from a single Sersic profile 
may be due, at least partially, to this extended envelope. 
This suggests that an additional model component may be 
required for them. We will re-visit two-component models 
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Figure 8. Distribution of the Sersic index n (upper left), effective 
radius Re (upper right), \ogRFFic (lower left) and logx? (lower 
right) from single Sersic fits for the BCGs divided by morphology. 
The red solid line corresponds to cD galaxies, the green dashed 
line to ellipticals, and the blue dotted line to spirals and SOs. 
The (j value in each panel indicates the significance (confidence 
level) of the observed differences between the cD and elliptical 
BCG parameter distributions. These are derived from two-sample 
Kolmogorov-Smirnov tests. 


elliptical BCGs than xS- 


Figure 6. The top four panels show the distribution of \ogRFF 
(left) and logx^ (right) for single Sersic (open histograms) and 
Sersic-hExponential (solid histogram) fits. The two uppermost 
panels correspond to the 53 1C BCGs, while the middle pan¬ 
els correspond to the 25 2C BCGs. The two bottom panels show 
the difference in RFF and between one-component and two- 
component models for both sets of BCGs. Clearly, the RFF 
and xi distributions of one- and two-component fits are vir¬ 
tually indistinguishable for 1C BCGs. However, RFF and x? 
tend to be significantly smaller for the two-component fits of 2C 
BCGs. Typical values for good fits are logiiFF~ —I.TI^q'Jq, and 
logxj — 0*042 ^q'q 25 (median +/— the 1st and 3rd quartiles of 
the parameters). Both RFF and x^ ^re sensitive to the magni¬ 
tude of the residuals, but RFF is appears to be significantly more 
sensitive. 


in Section EM The clear difference in RFF suggests that 
RFF could be another good discriminator to separate cD 
and elliptical BCGs. 

Similar conclusions can be reached from the the distri¬ 
butions of xi shown in the lower right panel of Figure [S] 
albeit less clearly. This is not surprising since, as shown in 
Section EM both RFF and xt measure the strength of the 
residuals, but xi is significantly less sensitive. Therefore, 
RFF is expected to be more efficient for separating cD and 


These results show a clear link between the visual morpholo¬ 
gies of BCGs and their structural properties. Although cD 
galaxies tend to have similar shapes to elliptical BCGs, they 
usually have larger sizes and their structures generally de¬ 
viate more from single Sersic profiles. In contrast, elliptical 
BCGs tend to be smaller, and their light profiles are sta¬ 
tistically more consistent with single Sersic models. These 
structural differences, especially in Re and RFF, could there¬ 
fore provide quantitative ways to separate elliptical and cD 
BCGs without relying on visual inspection. We will explore 
these issues in Section [5] 

4.3 Sersic+Exponential Models 

The RFF distributions shown in Section r4.2l indicate that el¬ 
liptical BCGs are statistically better fitted by a single Sersic 
model than cDs. Since a distinctive feature of cD galaxies is 
their extended luminous halo, two-component models may 
be more appropriate to describe accurately the light dis¬ 
tributions of cD BCGs. Following ISeigar et al.l ll200i1i and 
iDonzelli et al.l (l201lli . we explore here how a model consist¬ 
ing of an inner Sersic profile and an outer exponential enve¬ 
lope performs when fitting BCG images. The fitting process 
was described in detail in Section El 

As shown in Section El both RFF and xi can pro¬ 
vide quantitative information to assess whether BCGs are 
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Figure 9. Comparison of the residuals between single Sersic 
and Sersic+Exponential models. The left panel shows the frac¬ 
tional differences in RFF obtained with two-component and one- 
component fits for cD (red solid line), elliptical (green dashed 
line), and disk (blue dotted line) BCGs. The right panel shows 
the corresponding fractional differences for xS- 


better fitted by a two-component model than by a one- 
component model, at least in very clear cases. Figure [9] 
shows a comparison of these parameters obtained for sin¬ 
gle Sersic and Sersic-I-Exponential models. In the left panel 
we show a histogram of the fractional differences in the 
RFF values {RFFi^ - RFF 2 c)/RFFi^ for all three BCG 
types. The right panel shows the corresponding xi frac¬ 
tional differences (xfic — X?, 2 c)/x?,ic- ^ is clear that for 
disk BCGs, the Sersic-I-Exponential model does a better job. 
This is not surprising since spiral and lenticular galaxies con¬ 
tain clearly distinct bulges and disks. For elliptical BCGs 
the improvement in RFF and xS for two-component mod¬ 
els is generally quite small, as expected: elliptical galaxies 
are known to be reasonably well fitted by Sersic models, so 
the extra component does not improve the residuals signif¬ 
icantly. Perhaps surprisingly, the improvement is also only 
marginally better for cDs: the typical fractional differences 
for cD galaxies are {RFFic — RFF 2 c)/RFFic = 0.1llg;Jg 
and (xfic - xl, 2 c)lxl,ic = 0.035l!5;2|^ (median +/- 1st 
and 3rd quartiles of the parameter distributions). 

Since the distributions shown in Figure [9] for ellipti¬ 
cals and cDs are statistically indistinguishable, there is no 
clear separation that could be used to distinguish ellipti¬ 
cal and cD BCGs by comparing one-component and two- 
component fits. Moreover, on average, Sersic-I-Exponential 
model does not fit the profile of cD BCGs clearly better than 
single Sersic model. The reason is that for cD BCGs the 
values of RFF and x^ are generally not dominated by the 
presence or absence of a second exponential model compo¬ 
nent but by other structures present in the residual images, 
such as double cores. Since there is no clear improvement in 
the Sersic-I-Exponential model, the model with the smallest 
number of parameters (i.e., single Sersic model) will be pre¬ 
ferred for simplicity. The following discussions are based on 
the results from the single Sersic fits. 


4.4 Summary of Section [4] 

In this section we have analysed the differences in the 
structural properties of BGGs as a function of mor¬ 
phology. These structural parameters have been de¬ 


rived from one-component (Sersic) and two-component 
(Sersic-l-exponential) model fits. Disk BGGs (a small mi¬ 
nority) have smaller Sersic indices (n) than elliptical and 
cD BGGs, as expected. They also have different, generally 
broader, distributions of RFF and Xv- Elliptical and cD 
BGGs have similar n values, but cDs tend to have larger 
values of Re, RFF and Xv- These differences do not depend 
strongly on whether we use one- or two-component models. 

The observed structural differences could provide quan¬ 
titative ways to separate elliptical and cD BGGs without 
relying on visual inspection. We explore these in section [5] 
Furthermore, the differences we have found in the structural 
parameters suggest that the formation histories of elliptical 
and cD BGGs may be different. For instance, gas-rich major 
mergers and other dissipative processes may be responsible 
for building the inner (Sersic-like) component, while dissipa¬ 
tionless minor mergers may contribute to the build-up of the 
outer extended envelope and to the growth of galaxy sizes 
(e.g., lOser et ^ 120101 : 1 Johansson et al.l l2012l : iHuang et al.l 
I2OI3I) . We will explore in a subsequent paper (Zhao et ah, 
in preparation) whether the morphological and structural 
properties of BGGs are linked to other intrinsic BGG prop¬ 
erties such as their stellar mass, and/or to the properties 
of their environment. These links will provide more clues to 
the formation history of cDs/BGGs. 


5 SEPARATING ELLIPTICAL AND CD BCGS 

The results of Section suggest that we may be able to 
use the different distributions of cD and non-cD BGGs on 
the logRe-iogRFFie plane to separate them in an objec¬ 
tive, quantitative and automatic way. Figure [TU] shows that 
cDs are clearly segregated from other BGGs in this two- 
dimensional parameter space. We attempt to find a robust, 
well-defined way to separate, statistically, cD and non-cD 
BGGs using the information provided by this diagram. In 
other words, we suppose to find an “optimal border” that 
can separate them. 


5.1 Method Description and the Optimal Border 

Ideally, any process that selects cD galaxies from a sam¬ 
ple of BCGs needs to have high completeness (i.e., select as 
many of the cDs present in the sample as possible), while 
avoiding contamination from non-cDs (i.e., maximising the 
purity of the sample). These two requirements compete with 
each other, and increasing completeness often results in a de¬ 
crease in sample purity, and vice-versa. We need therefore to 
find the best compromise between these competing require¬ 
ments. In general, the optimal solution will depend on the 
specific intent for the selected sample, and therefore on the 
decision of how much weight to give to completeness and to 
purity. It is useful to define a measurement on the quality of 
the selection method that combines both requirements in a 
well-defined way. The optimal solution will then be obtained 
by maximising this quality param eter. 

Following IHovos et al.l ll201ll) the sensitivity, which is 
often known as completeness in astronomy, is defined as: 

_ #TruePositives , , 

/(:TruePositives -|- /i^FalseNegatives 
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Similarly, we define specificity as: 

_ #TrueNegatives , > 

^ ^TrueNegatives + ^FalsePositives 

A “True Positive” is an object retrieved by the selection 
process with the required properties (i.e., a cD galaxy that 
is correctly selected as such). A “False Negative” is an item 
that is not retrieved by the selection process but does present 
the needed properties (a cD galaxy that is not selected). A 
“True Negative” is an item that is rightfully rejected by the 
selection process since it does not have the required prop¬ 
erties (for instance, an elliptical galaxy that is not selected 
as a cD). A “False Positive” is an item that is incorrectly 
picked up by the selection process, but does not have the 
properties of interest (for example, an elliptical galaxy that 
is wrongly selected as a cD). 

Sensitivity and specificity can be combined into a single 
number, known as the F-score llvan Riisbergen|[l979l f . which 
provides a single measure on the quality of the selection 
process. The F-score is just a weighted harmonic average of 
r and p, 


^ (l-f/3^) xpxf 
^ 13^ xp + r 


( 6 ) 


where /3 is a control parameter that regulates the relative 
importance of completeness with respect to specificity. This 
is a user-supplied value that depends on the particular goals 
of the study. We will explore later how the choice of fi af¬ 
fects our selecting results. At this stage, a value of /? = 1.25 
is used, which can be thought of as weighing completeness 
more than the lack of contamination. For our BCG samples, 
the F-score is used to grade the performance of the diag¬ 
nostics we use when separating cD galaxies from the parent 
population. 

The selection process that we will apply to the parent 
population of BCGs in order to select cD galaxies will be de¬ 
fined by a “border” in the logFe-logFFFic plane (see Fig¬ 
ure [iQl). This border will be represented by a second-order 
polynomial in the horizontal coordinate. Higher-order poly¬ 
nomials (or more conplex functions) could be used, but the 
additional complexity is not required here. In our specific 
problem, the cD galaxies play the role of the “items pre¬ 
senting the required properties” discussed above, and the 
parent population is the complete sample of BCGs. 

According to the dehnition of sensitivity and specificity, 
the BGGs in the parent sample are classified into four cate¬ 
gories by their position relative to the border. In the logFe- 
logFFFic plane, cD galaxies dominate the region of large 
Fe and FFFic. We therefore define this region as the “cD 
side”. Thus 


• cD galaxies that fall on the cD side of the border are 
True Positives. 

• cD galaxies that do not fall on the cD side of the border 
are called False Negatives. 

• elliptical and disk (spiral and SO) BCGs that fall on the 
cD side are regarded as False Positives. 

• elliptical and disk (spiral and SO) BCGs that do not fall 
on the cD side of the border are True Negatives. 

The optimal border is found by maximising the F-score 
value. Following the method described in IHovos et al. 
ll201lll . we use the Amoeba algorithm (IPress fc Spergell 


3.5r 


3.0 


2.5 


F=0.69, r=0.75, p=0.61 
y=0.57+0.08x+0.31x^ 
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Figure 10. logi?e vs. log RFFic for the BCGs in our sample. We 
use this diagram to find the optimal border to separate cD from 
non-cD BCGs. The symbols are the same as in Figure[4] The black 
dotted line is the “first guess” for the border. The blue solid curve 
is the optimal border determined when we consider all cD BCGs 
(cD, cD/E and cD/SO) as cD galaxies. The blue dashed curve 
is the optimal border determined when we consider only pure 
cD and pure elliptical BCGs (excluding all cD/E, cD/SO, E/cD, 
E/SO, spiral and SO BCGs). The legend shows the maximum F- 
score for the optimal borders and the corresponding completeness 
r and specificity p. The equations defining the optimal borders are 
also shown. The error bar shows the mean error of each parameter. 
We used /3 = 1.25 in this case. 


llQSSh to carry out this maximization and find the polyno¬ 
mial defining the border. 

It is clear from Figure^Jthat the selected galaxy sam¬ 
ple on the cD side of the optimal border will not contain only 
cD galaxies, and a degree of contamination will be present. 
We define contamination dHovos et al.ll201ll 'l as: 

^ #non-cDs tested as positive 

G ~ - 7 , —7^-:—:-• ( I ) 

#afl positives 

The numerator are the non-cD BCGs which are on the cD 
side of the optimal border. The denominator of this fraction 
includes both cD galaxies and non-cD BCGs on the cD side. 

Figure [Til] shows the log Fe-log FFFic plane for the 
BCGs in our sample. The Amoeba algorithm requires a first 
guess for the border, shown by the black horizontal dotted 
line. The optimal border determined by the algorithm does 
not depend on the exact initial guess. The blue solid curve 
is the optimal border determined when we consider all cD 
galaxies (cD, cD/E and cD/SO) as cD galaxies. This bor¬ 
der, computed using /3 = 1.25, has F-score = 0.69. 75% of 
all the cD galaxies are above the border (r = 0.75), and 
thus selected from the parent sample. The remaining 25% 
are mixed with the elliptical and disk BCGs in the region 
below the border. This selection therefore yields 75% com¬ 
pleteness. The galaxy sample above the border contains 311 
cD galaxies and 79 non-cD BCGs resulting in a ~ 20% con¬ 
tamination in the selected cD samples. In the region below 
the border there are 103 cD galaxies and 122 and non-cD 
BCGs. Thus, the non-cD BGC sample has a contamination 
of 46% from cD galaxies. This indicates that this technique 
is more effective (cleaner) at selecting cD galaxies than at 
selecting non-cD BCGs. 
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Figure 11. Two-step process to select cD BCGs. Symbols and 
legend are the same as in Figure fTOl Disk (spiral and SO) BCGs 
are separated from non-disk BCGs (cDs and ellipticals) first using 
the optimal border shown as the blue dashed curve. cD galaxies 
are then selected using the optimal border shown as the blue solid 
curve. See text for details. 



a completeness of 77% (305 cDs are selected), and a con¬ 
tamination of only 14%. Compared to the single-step cD 
selection (311 cDs were selected with 20% contamination), 
the two-step process clearly selects a very similar number 
of cDs but with better purity. The decision on whether the 
increase in purity is worth the additional complexity is left 
to the reader. In the reminder of this paper we will use the 
single-step selection process for simplicity. 

The automatic techniques we have developed can be 
applied to any BCG sample, but the optimal border needs 
to be adapted and calibrated using the imaging data from 
which the parent sample was derived. The calibration can 
be performed using a sub-sample of visually-classified BCGs, 
and then automatically applied to the complete sample using 
the structural parameters determined from standard single- 
Sersic fits. 

A /3 value needs to be chosen depending on whether we 
are more interested in the completeness of the cD sample 
or in its purity, but we suggest that P = 1.25 represents a 
reasonable compromise (see section [Oil . Furthermore, it is 
important to remember that this method works better at 
selecting a sample of cD galaxies rather than a sample of 
non-cDs. 


Note that if we consider a “cleaner” sample that con¬ 
tains only pure cD and pure elliptical BCGs (excluding all 
cD/E, cD/SO, E/cD, E/SO, spiral and SO BCGs), the opti¬ 
mal border (blue dashed curve in Figure[TD]) does not change 
significantly, but the quality of the selection as determined 
by the E-score value, the completeness r and the specificity 
p improves. This is not surprising: the identification of BGGs 
as pure cDs/Es (as opposed to the “dubious” ones) depends 
on more secure morphological characteristics which should 
be linked more clearly to the structural parameters. How¬ 
ever, considering only this cleaner sample is not a realistic 
scenario since in practical cases we would like to start from 
a full sample of BCGs and find which ones are cDs. Never¬ 
theless, it is reassuring that the border we determine does 
not depend very strongly on the exact training set used. 

On the selected cD side, spiral BGGs are an important 
source of contamination. However, since most of them ap¬ 
pear in the large RFFic region, it would be possible to go 
a step further to implement a simple further refinement in 
our method to separate spirals from the selected cDs: very 
few cD galaxies have log RFFic larger than ~ —1.1. This 
would significantly improve the purity of the cD sample at 
very little cost in terms of its completeness. 

Moreover, it is clear from Figure [in] that all disk BGGs 
(spirals and SOs) contribute significantly to the contamina¬ 
tion of either the cD or the elliptical samples separated 
by the best border. However, we can use the fact that 
disk BGGs distribute over a distinct area on the log Re- 
log RFFic plane to apply a two-step process to exclude them 
from our cD selection. First, the disk BGGs can be separated 
from the elliptical and cD BCGs, and then the cD BGGs 
can be selected out of the rest BCG sample. Figure [11] illus¬ 
trates the results of this two-step selection. The blue dashed 
curve is the optimal border determined in the the first step. 
By excluding disk BGGs using this border, a very complete 
(r = 0.93) and pure (p = 0.87) non-disk BGG sample is 
built. The cDs can then be separated from the ellipticals 
using the optimal border shown by the blue solid curve with 


5.2 Distance to the Optimal Border 

It is informative to explore the distribution of the points 
in the log i?e-log EEEic plane (Figure fTUl) in terms of their 
minimum (perpendicular) distance to the optimal border. 
We define the distance from each point to the optimal border 
as 

D = /(-AlogEEEr.y^/Alo^y^ 

V V O-logflFFie / V <^108 fle / 

where Alog/ZEEic is the difference in log RFFic between 
the data point and the optimal border, and uiogijFFio is the 
dispersion in log RFFic computed for all the points. A log i?e 
and (Tiog have a similar meaning but for log Re. Note that, 
because the units of the x and y axes are different, the dis¬ 
tance is measured in units of the scatter of each parameter. 
For each point, the minimum distance Dmin can be then 
determined. Figure [T^ shows the distribution of these mini¬ 
mum distances for the different morphologies. As expected, 
the vast majority (> 80%) of the cDs show positive distances 
(they are above the optimal border line) while most of the 
ellipticals have negative ones. Under 20% of the cDs spill 
over to the negative region, severely contaminating the non- 
cD sample, while a few ellipticals weakly contaminate the 
cD region. The measurement errors in log Re (~ 0.16) and 
log RFFie (~ 0.13) result in distance errors on the order of 
0.7 in this metric. This contributes to the cDs’ “spillover”, 
but does not completely explain it. Reducing the measure¬ 
ment errors would certainly improve the performance of our 
method, but it would never make it perfect. 

Interestingly, the spiral and SO BGGs are quite well sep¬ 
arated: the former show mostly positive distances while the 
later have mostly negative ones. This is mainly due to spi¬ 
rals having generally larger RF Fie values because the spiral 
arms and star-forming regions are not included in the Sersic 
models, while the SOs are smoother. This clear separation 
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Figure 12. Distribution of the minimum distances to the opti¬ 
mal border shown in Figure [TOl for the cD and elliptical BCGs 
(top panel) and the spiral and SO BCGs (bottom panel). Positive 
and negative distances correspond to points above and below the 
optimal border line respectively. 


provides a possible way to separate spiral and SO galaxies, 
but this needs to be further tested with large disk samples. 

Another interesting result is that BCGs classified as 
pure and uncertain cDs (e.g., cD/E) have very different min¬ 
imum distance distributions (Figure [13] top panel). About 
half of the cD/E BCGs have negative distances (i.e., are on 
the wrong side of the border), but only ~ 20% of the pure 
cDs do. Most of the spillover of the pure cDs into the nega¬ 
tive region, however, can be explained by the measurement 
errors. It should be noticed that the difficulties inherit in 
the visual morphological classification are directly reflected 
in the structural parameters: when the visual classifier is cer¬ 
tain that a BGG is a cD, its structural parameters almost 
always confirm it, while in uncertain cases (e.g., cD/E) the 
structural parameters reflect this uncertainty. Similar con¬ 
clusions can also be obtained from the pure elliptical BCGs 
and uncertain ones (e.g., E/cD), as shown in the bottom 
panel of Figure [T^ 

This analysis confirms the visual impression in terms of 
the BCG structure that there is a continuous distribution in 
the properties of the BGG extended envelopes, ranging from 
undetected (pure E class) to clearly detected (pure cD class), 
with the intermediate classes (E/cD and cD/E) showing in¬ 
creasing degrees of envelope presence. This continuous dis¬ 


Figure 13. Distribution of the minimum distances to the opti¬ 
mal border shown in Figure [TOl for the pure cD BCGs and cD/E 
BCGs (top panel). The bottom panel shows the corresponding 
histograms for pure E BCGs and E/cD BCGs. 

tribution in envelope detectability is reflected quantitatively 
in the structural parameters of the BGGs, by the minimum 
distance to the optimal border providing some indication of 
the relative importance of the envelope. 

5.3 Effect of the jS Parameter 

In the E-score definition, the /3 parameter is used to ap¬ 
portion weight to the completeness and the specificity. For 
larger values of P the completeness is given a larger weight 
than the lack of contamination. Conversely, smaller values 
of P prioritise lack of contamination above completeness. 
To test how changing P affects the results of the selection 
process, we repeat the exercise carried out in Section 15.11 
but using P = 2.0 and P = 0.5 in the determination of the 
optimal border. 

Figure [TT] shows the optimal border for p = 2.0 (upper 
panel) and P = 0.5 (lower panel). It is clear that the P 
parameter has a decisive impact on the selection of potential 
cD galaxies. As shown in the upper panel, when compared 
to the P = 1.25 results, 11% more galaxies are correctly 
identified as cDs, significantly increasing the completeness. 
The price paid is that the specificity goes down from 61% 
to 46% since more non-cD BCGs are included. Gonversely, 
in the lower panel (/3 = 0.5) the selected cD sample is purer 

























































The Link Between Morphology and Structure of BCGs 15 



log RFFi, 



log RFFi, 


Figure 14. Illustration of the effect of /3 on the optimal border. 
The symbols, lines and legends have the same meaning as in Fig¬ 
ure llUl but we use /3 = 2.0 for the upper panel and /3 = 0.5 for the 
lower panel. With /3 = 2.0 we give more weight to the complete¬ 
ness than to the lack of contamination. When using 0 = 0.5, the 
lack of contamination is given more importance than achieving 
higher completeness. The choice on 0 depends on the aims of the 
specific research. 


{p = 0.85), but at the expense of completeness, with 20% 
fewer cD galaxies selected when compared with the 0 = 1.25 
result. 

With 0 = 2.0, the contamination of the cD sample by 
non-cDs is 23%, while the contamination of the non-cD sam¬ 
ple by cDs is 39%. With 0 = 0.5, the corresponding values 
are 12% and 52% respectively. Therefore, for any value of 
0 this selecting technique is cleaner and more effective at 
selecting cD galaxies than at selecting non-cD BCGs. 

As before, if we consider a cleaner sample that contains 
only pure cD and pure elliptical BCGs, the optimal border 
(blue dashed curve) does not change significantly, but the 
T-score value, the completeness r and the specificity p im¬ 
prove. However, we have argued that this does not represent 
a realistic scenario. 

We conclude that 0 = 1.25 represents a good compro¬ 
mise, as its optimal border picks up a cD galaxy sample 
reasonably complete, and with relatively small contamina¬ 
tion. However, no single value of 0 can be considered to be 


“correct” and needs to be set according to the scientific goals 
of the study. 


6 CONCLUSIONS 

In this paper we have analysed a well-defined sample of 
625 low-redshift Brightest Cluster Galaxies published in 
Ivon der Linden et al. ll2007l ) with the aim of linking their 
morphologies to their structural properties. We morpholog¬ 
ically classified the BCGs using SDSS r-band images and 
found that over half of them (~ 57%) are pure cD galaxies 
and pure elliptical BCGs constitute ~ 13% of the sample. 
The intermediate classes (mostly cD/E or E/cD) account 
for ~ 21%. It suggests a continuous distribution in the prop¬ 
erties of the BGG extended envelopes, ranging from unde¬ 
tected (pure E class) to clearly detected (pure cD class), with 
the intermediate classes (E/cD and cD/E) showing increas¬ 
ing degrees of envelope presence. We found this continuous 
distribution in envelope detectability is reflected quantita¬ 
tively in the structural parameters of the BCGs. There is 
also a minority of BCGs that are neither cD nor ellipti¬ 
cal. About 7% are disk galaxies (spirals and SOs, in similar 
proportions) and the rest (~ 2%) are in merging (see ap¬ 
pendix A). 

In order to link the morphologies of the BCGs to their 
structural properties, we have fitted the BCGs light distri¬ 
butions with the SDSS r-band images using one-component 
(Sersic) and two-component (Sersic-1-Exponential) models. 
We first characterised how well the models fit the target 
BCG by using two quantitative diagnostics. One diagnostic 
is the residual flux fraction (RFF), which measures the frac¬ 
tion of the galaxy flux presenting in the residual images after 
subtracting the models. The other diagnostic is the reduced 
Xu- We concluded that generally it is very difficult to find 
a robust diagnostic to decide, in a statistic way, whether a 
one-component or a two-component model is preferred for 
BCGs, especially for cD galaxies. Since there is no evident 
improvement by using two-component model fits, our other 
conclusions rely on the one-component Sersic fits. 

From simple one-component Sersic profile fits, we have 
found a clear link between the BCGs morphologies and their 
structures, and claimed that a combination of the best-fit pa¬ 
rameters can be used to separate cD galaxies from non-cD 
BCGs. In particular, cDs and non-cDs show very different 
distributions in the Re-RFFic plane, where Re is the ef¬ 
fective radius and RFFie is the residual flux fraction, both 
determined from Sersic fits. cDs have, generally, larger Re 
and RFFic values than ellipticals. Therefore we found, in 
a statistically robust way, a boundary to separate cD and 
non-cD BCGs in this parameter space. BCGs with cD mor¬ 
phology can be selected with reasonably high completeness 
(~ 75%) and low contamination (~ 20%). 

This automatic and objective technique can be applied 
to any current or future BCG samples which have good qual¬ 
ity images. The method needs to be adapted and calibrated 
using the imaging data from which the parent sample was 
derived. Once calibrated with a representative sub-sample 
of visually-classified BCGs, this technique can be applied to 
the complete sample using the structural parameters deter¬ 
mined from standard single-Sersic fits. 

In a subsequent paper (Zhao et ah, in preparation) we 
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will explore how the morphological and structural properties 
of BCGs are linked to other intrinsic BCG properties such as 
their stellar mass, and/or to the properties of their environ¬ 
ments. These links will provide more clues to the formation 
history of cDs/BCGs. 
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APPENDIX A: DATA TABLE 

Table A1 contains the main properties of the BCGs dis¬ 
cussed in this paper. The full table is published electroni¬ 
cally. 


ID2 

(1) 

ID3 

(2) 

RA 
deg (3) 

DEC 
deg (4) 

Z 

(5) 

•^cl 

kms“l (6) 

logRe.lc 
kpc (7) 

nic 

(8) 

RFFic 

(9) 

X?c 

(10) 

Type 

(11) 

Comments 

(12) 

1011 

1013 

227.107346 

-0.266291 

0.091 

748 

1.527 

5.38 

0.08190 

1.752 

cD 

Clear halo; perhaps interacting 

1023 

1025 

153.409478 

-0.925413 

0.045 

790 

1.908 

6.25 

0.05052 

1.374 

cD 

Clear halo; interacting with fainter galaxies 

1064 

1075 

153.437067 

-0.120224 

0.094 

875 

1.312 

4.49 

0.02648 

1.086 

E/cD 


- 

1027 

191.926938 

-0.137254 

0.088 

1020 

1.063 

4.42 

0.06594 

1.903 

E 

Interacting/merging with bright early-type 

- 

1389 

202.337884 

0.749685 

0.080 

853 

1.044 

6.02 

0.01990 

1.087 

E/cD 

Faint/small halo 

2040 

2050 

17.513187 

13.978117 

0.059 

759 

2.408 

9.77 

0.04122 

1.224 

cD 

Several bright-ish companions 

1052 

1058 

195.719058 

-2.516350 

0.083 

749 

1.627 

4.89 

0.04694 

1.455 

cD 

Multiple merger 

1034 

1036 

192.308670 

-1.687394 

0.085 

771 

0.977 

4.86 

0.02102 

1.115 

E 


1041 

1044 

194.672887 

-1.761463 

0.084 

771 

2.318 

5.64 

0.05023 

1.280 

cD 

Very large, elongated halo; some faint companions 

- 

1126 

192.516071 

-1.540383 

0.084 

878 

2.039 

9.12 

0.04520 

1.348 

cD 

Interacting with faint companions 

3002 

3004 

258.120056 

64.060761 

0.080 

1156 

1.667 

4.81 

0.02561 

0.991 

cD 


3096 

3283 

135.322540 

58.279747 

0.098 

756 

1.866 

6.96 

0.05535 

1.144 

cD 

Merging with bright companion 

1045 

1048 

205.540176 

2.227213 

0.077 

828 

0.883 

2.52 

0.10689 

11.280 

E/cD 

Multiple merger 

1003 

1004 

184.421356 

3.655806 

0.077 

966 

1.753 

4.75 

0.05233 

1.225 

cD/E 

Interacting/merging with early-type 

- 

1456 

173.336242 

2.199054 

0.099 

746 

1.696 

8.09 

0.02573 

1.128 

cD 


1053 

1061 

228.220703 

4.514004 

0.038 

789 

0.875 

7.54 

0.01749 

1.074 

cD 


2163 

2074 

314.975446 

-7.260758 

0.079 

765 

1.231 

8.03 

0.04481 

1.315 

E/cD 


2002 

2002 

358.557007 

-10.419200 

0.076 

812 

2.660 

11.12 

0.03832 

1.201 

cD 

Many faint and bright-ish companions 

2006 

2013 

10.460272 

-9.303146 

0.056 

903 

1.433 

1.62 

0.04140 

1.477 

cD 

Several faint companions 

1355 

1460 

175.554108 

5.251709 

0.097 

1074 

0.952 

5.30 

0.01557 

1.052 

cD 

Interacting with faint galaxy; faint but clear halo 

1058 

1069 

184.718166 

5.245665 

0.076 

721 

1.988 

7.98 

0.04144 

1.251 

cD 

Interacting with faint galaxies 

1002 

1002 

159.777581 

5.209775 

0.069 

800 

1.740 

8.40 

0.03838 

1.321 

cD/E 

Clear halo 

- 

1276 

183.271286 

5.689677 

0.081 

729 

0.995 

5.30 

0.02142 

1.151 

E 


1039 

1042 

228.808792 

4.386210 

0.098 

857 

1.800 

8.77 

0.04365 

1.205 

E/cD 

Some halo? faint companions 

- 

3332 

124.471428 

40.726395 

0.063 

802 

1.463 

6.40 

0.08125 

2.639 

SBO 


3011 

3028 

204.034694 

59.206401 

0.070 

872 

2.120 

7.86 

0.08172 

1.556 

cD 

Several faint companions 

1001 

1001 

208.276672 

5.149740 

0.079 

746 

1.820 

7.85 

0.02720 

1.128 

E/cD 


3004 

3012 

255.677078 

34.060024 

0.099 

1127 

1.717 

3.54 

0.08433 

1.949 

cD 

Late merger? 

- 

3094 

254.933115 

32.615319 

0.098 

875 

1.291 

3.50 

0.02878 

1.069 

cD 

Very faint companions 

- 

1066 

202.795126 

-1.730259 

0.085 

814 

1.942 

9.09 

0.03653 

1.161 

E/cD 

Interacting/merging with bright galaxy and fainter one 

- 

2214 

321.599487 

10.777511 

0.095 

741 

0.818 

3.98 

0.02260 

1.199 

E 


2096 

2109 

359.836166 

14.670211 

0.093 

786 

1.161 

6.56 

0.03572 

1.242 

cD/E 


2085 

2085 

334.197449 

-9.724778 

0.094 

806 

0.779 

3.43 

0.02861 

1.348 

cD 


2027 

2035 

4.177309 

-0.445436 

0.065 

1084 

1.436 

8.89 

0.02417 

1.168 

cD 

Several companions 

- 

3084 

118.360820 

29.359459 

0.061 

781 

1.584 

3.95 

0.06632 

1.382 

cD 

Several faint and bright companions 

- 

3347 

119.679733 

30.773809 

0.076 

902 

1.354 

6.04 

0.01470 

1.019 

E/cD 


- 

1283 

125.745443 

4.299105 

0.095 

754 

2.747 

10.47 

0.04483 

1.094 

cD 

Several faint-ish companions 

- 

1039 

186.878093 

8.824560 

0.090 

846 

1.962 

6.94 

0.06100 

1.965 

cD 

Clear halo, bright companion (dumbbell galaxy) 


Table Al. Properties of the BCG sample. Columns (1) and (2) provide sa’la-xy iden tifications, where ID2 is the SDSS-C4 number <SDSS-C4 NNNN> and IDS is the SDSS C4_2003 
number, <SDSS-C4-DR3 NNNN>, as given in Simbad (Ivon der Linden et alJuOOTh . Columns (3) and (4) give the right ascension and declination in degrees. Column (5) gives the 
redshift and column (6) the velocity dispersion of the cluster. Columns (7), (8), (9) and (10) contain the effective radius, Sersic index, residual flux fraction and reduced derived from 
the single Sersic fits (see text for details). Column (11) gives the visual morphological classification of the BCGs. Column (12) contains some comments from the classifier. 
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